SUMS OF LAPLACE EIGENVALUES ROTATIONALLY 
SYMMETRIC MAXIMIZERS IN THE PLANE 



R. S. LAUGESEN AND B. A. SIUDEJA 

Abstract. The sum of the first n > 1 eigenvalues of the Laplacian 
is shown to be maximal among triangles for the equilateral triangle, 
maximal among parallelograms for the square, and maximal among el- 
lipses for the disk, provided the ratio (area) ^/(moment of inertia) for 
the domain is fixed. This result holds for both Dirichlet and Neumann 
eigenvalues, and similar conclusions are derived for Robin boundary con- 
ditions and Schrodinger eigenvalues of potentials that grow at infinity. 
A key ingredient in the method is the tight frame property of the roots 
of unity. 

For general convex plane domains, the disk is conjectured to maxi- 
mize sums of Neumann eigenvalues. 



1. Introduction 

Eigenvalues of the Laplacian represent frequencies in wave motion, rates 
of decay in diffusion, and energy levels in quantum mechanics. Eigenvalues 
are challenging to understand: they are known in closed form on only a 
handful of domains. This difficulty has motivated considerable work on 
estimating eigenvalues in terms of simpler, geometric quantities such as area 
and perimeter. 

We will obtain a sharp bound on the sum of the first n > 1 eigenvalues 
of linear images of rotationally symmetric domains. Our methods apply 
equally well to Dirichlet, Robin, and Neumann boundary conditions. 

Write Ai,A2,A3,... for the Dirichlet eigenvalues of the Laplacian on a 
plane domain of area A and moment of inertia /. We will prove that for 
each n > 1, the normalized, scale invariant eigenvalue sum 

(Ai + --- + A„)^ (LI) 

is maximal among triangular domains for the equilateral triangle. Among 
parallelograms the maximizer is the square, and the disk is the maximizer 
among ellipses. The only case known previously was the fundamental tone, 
n = 1, due to Polya [44]. 
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An analogous result will be shown for the sum of Neumann eigenvalues, 

^3 

(//I + ••• + //„) — , (1.2) 

and then for Robin and Schodinger eigenvalues too. These latter results are 
new even for n = 1. See Section [3l 

Our work suggests conjectures for general convex domains. Is the Dirichlet 
eigenvalue sum (jl.ip maximal for the disk? Not when n = 1, curiously, 
because any rectangle or equilateral triangle gives a larger value for the 
fundamental tone. We conjecture that those domains maximize XiA^/I. 
For the Neumann eigenvalue sum (II. 2p it does seem plausible to conjecture 
maximality for the disk, as we discuss in Section HI 

Central to the paper is a new technique we call the 

"Method of Rotations and Tight Frames" . 

The idea is to linearly transplant the eigenfunctions of the extremal domain 
and then average with respect to allowable rotations of that domain. This 
averaging of the Rayleigh quotient is accomplished using a "tight frame" or 
Parseval identity for the root-of-unity vectors. The Hilbert-Schmidt norm 
of the linear transformation arises naturally in such averaging, and is then 
represented in terms of the moment of inertia. 

Intuition. The eigenvalue sum (jl.ip can be written as a product of two 
scale invariant factors, as 

{\l + --- + Xn)A-—. 

The first factor, (Ai + • • • + A^)^, normalizes by the area of the domain, and 
so can be thought of as a generalized "Faber-Krahn" term (although the 
Faber-Krahn theorem says XiA is minimal for the disk, not maximal). 

The second factor, A'^/I, is purely geometric and measures the "deviation 
from roundness" of the domain. This factor is small when the domain is 
elongated, and hence it balances the largeness of the first factor on such 
domains. 

The motivating intuition is that for a domain with characteristic length 
scales a and b, we have 

^ - 62 {abf - A3 ■ 

This rough calculation is exact for rectangles, up to constant factors. 

The first inequality in the calculation extends readily to higher dimen- 
sions, but the algebraic identity in the middle becomes more complicated. 
Thus in higher dimensions it seems that the moment of inertia should be 
evaluated instead on some kind of "reciprocal" domain that has length scales 
1/a and 1/6 and so on. Higher dimensional results of this nature will be de- 
veloped in a later paper |35j; the maximizing domains are regular tetrahedra, 
cubes, and other Platonic solids. 
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Related work. The major contribution of this paper is that it proves upper 
bounds that are geometricahy sharp, on eigenvalue sums of arbitrary length. 
We do not know any similar results in the literature. 

Some results of a different type are known, as we now describe. A bound 
due to Kroger [26] for Neumann eigenvalues says that (//i + • • • + Hn)A/n'^ < 
2tt. The inequality is asymptotically sharp for each domain, because /U„A ~ 
47rn by the Weyl asymptotics. But Kroger's bound is not geometrically 
sharp for fixed n, because there are no domains for which equality holds. 
Kroger's result should be viewed as a weak version of the Polya conjecture. 
That conjecture asserts that the Weyl asymptotic estimate is in fact a strict 
upper bound on each Neumann eigenvalue. It has been proved for tiling 
domains by Kellner [23], and up through the third eigenvalue for simply 
connected plane domains by Girouard, Nadirashvili and Polterovich [18j, 
but it remains open in general. 

Kroger also proved an upper bound on Dirichlet eigenvalue sums, involv- 
ing e-neighborhoods of the boundary [27]. This bound is again not geometri- 
cally sharp. Kroger's estimates were generalized to domains in homogeneous 
spaces by Strichartz [53]. 

Weak versions of Polya's conjectured lower bound for Dirichlet eigenvalues 
[46| are due to Berezin [8j and Li and Yau [37j . with later developments by 
Laptev [29] and others using Riesz means and "universal" inequalities, as 
surveyed by Ashbaugh [4j. Useful upper bounds on eigenvalue sums in terms 
of other eigenvalue sums have lately been obtained this way, by Harrell and 
Hermi [201 Corollary 3.1]. Note Polya's lower bound has been investigated 
also for eigenvalues under a constant magnetic field, by Frank, Loss and 
Weidl [12]. 

There is considerable literature on low eigenvalues of domains constrained 
by perimeter, in-radius, or conformal mapping radius, rather than moment 
of inertia. We summarize this literature in Section |8l 

Eigenvalues of triangular domains have been studied a lot, in recent years 
PiaiiaiMlllIllMllMllMllMlllI], and this paper extends the theory of 
their upper bounds. Lower bounds on Dirichlet eigenvalues of triangles are 
proved in a companion paper [34]: there the triangles are normalized by 
diameter (rather than area and moment of inertia) and equilateral triangles 
are shown to minimize (rather than maximize) the eigenvalue sums. 

For broad surveys of isoperimetric eigenvalue inequalities, one can consult 
the monographs of Bandle [7], Henrot [21], Kesavan [25| and Polya-Szego 
[38] , and the survey paper by Ashbaugh [3]. 

2. Assumptions and notation 

Eigenvalues. For a bounded plane domain D, we denote the Dirichlet 
eigenvalues of the Laplacian by Xj{D), the Robin eigenvalues by pj{D;a) 
where the constant cr > is the Robin parameter, and the Neumann eigen- 
values by iJ,j{D). In the Robin and Neumann cases we make the standing 
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assumption that the domain has Lipschitz boundary, so that the spectra are 
well defined. Denoting the eigenfunctions by Uj in each case, we have 

{—^Uj = XjUj in D j = pjUj in D J — Auj = fJ-jUj in D 

Uj = ondD \^ + (TUj = ondD \ ^ = on dD 

and 

< Ai < A2 < A3 < . . . < pi < P2 < P3 < • • • = ^1 < /i2 < ^3 < 

The Robin case reduces to Neumann when a = 0, and formally reduces to 
the Dirichlet case when o" = 00. 

The corresponding Rayleigh quotients are 

Dirichlet: R[u] = ' ' , — for u G H},{D), 

Jj-, dx 

Robin: R\u] = ^ — for u G H\D), 

dx 

Neumann: R[u] = ' , — for u € H\D). 

dx 

The Rayleigh— Poincare Principle [TJ p. 98] characterizes the sum of the 

first n > 1 eigenvalues as 

Ai + • • • + A„ 

= min + • • • + R[vn] : vi, . . . ,Vn € Hq{D) are pairwise orthogonal in L^(Z))} 

in the Dirichlet case, and similarly in the Robin and Neumann cases (using 
trial functions in instead of H^). 

Geometric quantities. Let 

A = area, 

/ = moment of inertia about the centroid. 

That is, 

I{D) = / \x-x\'^dx 
Jd 

where the centroid is x = ^^pjy Jjjxdx. 

Given a matrix M, write its Hilbert-Schmidt norm as 



where denotes the transposed matrix. 
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3. Sharp upper bounds on eigenvalue sums 

Dirichlet and Neumann eigenvalues. Our first result examines tlie ef- 
fect on eigenvalues of linearly transforming a rotationally symmetric domain, 
like in Figure [1] 

Theorem 3.1. If D has rotational symmetry of order greater than or equal 
to 3, then 

(Ai + • • • + A„)|^(^) < 1\\T-Yhs (Ai + • • • + A„)|^ (3.1) 

for each n > 1 and each invertible linear transformation T of the plane. The 
same inequality holds for the Neumann eigenvalues. 

Equality holds in (|3.1|) for the first Dirichlet eigenvalue (n = 1) if and 
only if either 

(i) T is a scalar multiple of an orthogonal matrix, or 

(ii) D is a square and T{D) is a rectangle (possibly with sets of capacity 
zero removed). 

Equality holds for the second Neumann eigenvalue (n = 2) if and only if T 
is a scalar multiple of an orthogonal matrix. 

The proof is in Section [6l Notice equality does hold in the theorem when 
T is a scalar multiple of an orthogonal matrix, because if T = rS" where S 
is orthogonal, then Aj(T(Z))) = r~'^\j{D) by rescaling and rotation, and 

1||T-1||2 _ „-2 

2ll^ Whs — r ■ 

The rotationally symmetric domain D in the theorem need not be convex, 
need not be a regular polygon, and need not have any axis of symmetry. For 
example, it could be shaped like a three-bladed propeller. 

Polya obtained the theorem for n = 1 (the Dirichlet fundamental tone), 
although with no equality statement. He stated this result in [44], and Polya 
and Schiffer proved it along with results for torsional rigidity and capacity 
in \4:7\ Chapter IV]. Our method differs subtly from theirs, as we explain 
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in Section [6l and this difference allows us to handle higher eigenvalue sums 
and Neumann eigenvalues too. 

To express the theorem more geometrically, we observe that the Hilbert- 
Schmidt norm of the transformation can be expressed in terms of mo- 
ment of inertia and area (Lemma l5.3p . Hence in particular: 

Corollary 3.2. Among triangles, the normalized Dirichlet eigenvalue sum 

(Ai + --- + A„)^ (3.2) 

is maximal for the equilateral triangle, for each n > 1. When n = 1, every 
maximizer is equilateral. 

Among parallelograms, the quantity (j3.2p is maximal for the square. When 
n = 1, every maximizer is a rectangle and every rectangle is a maximizer. 

Among ellipses, the quantity (j3.2|) is maximal for the disk. When n = 1, 
every maximizer is a disk. 

The normalized Neumann eigenvalue sum 

A^ 

(/i2 + ---+Mn) — 

is maximal among triangles for the equilateral triangle, among parallelo- 
grams for the square, and among ellipses for the disk, for each n> 2. When 
n = 2, every maximizer is an equilateral triangle, square, or disk, respec- 
tively. 

The Neumann case with ra = 1 is not interesting, because the first eigen- 
value equals for each domain. 

Remarks. 

1. The method extends to linear images of regular A^-gons for any N, 
but the most interesting cases are triangles and parallelograms (N = 3 and 
A'^ = 4), as considered in the corollary. 

2. For triangles, the moment of inertia can be calculated in terms of the 
side lengths h,l2,h as 

1=^(11 + 11 + 11)- (3.3) 

For a parallelogram with adjacent side lengths li,l2, the moment of inertia 
equals 

I=^{ll + ll)- (3.4) 

3. The eigenvalues and eigenfunctions of the extremal domains (the equi- 
lateral triangle, square and disk) are not used in our proofs. The eigenvalues 
are stated anyway in Appendix [Aj for reference. It is interesting to substi- 
tute them into Corollary 13.21 and obtain explicit estimates on eigenvalue 
sums. For example, for the Dirichlet fundamental tone (n = 1) one obtains 
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that 

127r for triangles, with equality for equilaterals, 

127r^ for parallelograms, with equality for rectangles, 

ill' 



2jQ^7r^ — ll.Svr^ for ellipses, with equality for disks. 



^3 

y«2— < 



(3.5) 

All three inequalities were obtained by Polya [S], |47l p. 308,328]. The 
first inequality, for triangles, was rediscovered with a different proof by Fre- 
itas |131 Theorem 1]. The second inequality was rediscovered for the special 
case of rhombi by Hooker and Protter [23, §5] and for all parallelograms by 
Hersch [22. formula (5)], again with different proofs. These authors stated 
their results in terms of side lengths, as 

(li +I2+ ^3)^73^' for triangles, 
~ 1 (^1 + ^D^^M^ parallelograms. 

These inequalities are equivalent to Polya's by formulas (j3.3|) and (j3.4|) for 
the moment of inertia. 

For the first nonzero Neumann eigenvalue we find from Corollary 13.21 and 
Appendix |A] that 

47r^ for triangles, with equality for equilaterals, 

67r^ for parallelograms, with equality for squares, 

2(j[ — G.Bvr^ for ellipses, with equality for disks. 

(3.7) 

These inequalities too can be stated in terms of side lengths. For stronger 
inequalities on fj.2 , see Section [HI 

For n = 3, the Corollary says {fi2+ fJ-s)^^ / 1 is maximal for the equilateral 
triangle. This result was proved recently by the authors using a different 
method with explicit trial functions |32[ Theorem 3.5]. 

4. Corollary 13.21 becomes false when applied to individual eigenvalues 
instead of eigenvalue sums. For example, X^A"^/! is not maximal for the 
square among rectangles: to the contrary, it is locally minimal. The under- 
lying reason is that a square has a double eigenvalue A2 = A3 that "splits" 
when the square is deformed into a rectangle of the same area; the second 
eigenvalue decreases and the third increases, while the moment of inertia 
varies only at second order. 

Robin eigenvalues. In the next theorem we normalize the Robin param- 
eter in terms of T~^, in order to obtain a scale invariant expression. 

Theorem 3.3. If D has rotational symmetry of order greater than or equal 
to 3, then 

< (piH + Pn)-j- 

<T\\T-mHs/V2,TiD) ^ a,D 

for each n > 1 and each invertible linear transformation T of the plane. 



A^ 

{Pl + ■ ■ ■ + Pn) — 
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Equality holds for the first Robin eigenvalue (n = 1) if and only if T is a 
scalar multiple of an orthogonal matrix. 

The subscript "o", D" on the right side of the inequahty specifies the do- 
main where the eigenvalues and geometric quantities are to be evaluated, 
and also the value of the Robin parameter to be used. The subscript on the 
left side of the inequality similarly specifies the domain T{D) and the Robin 
parameter a\\T~^\\HS / to be used there. 

Corollary 3.4. Fix the Robin parameter a > 0. Among all triangles of the 
same area, the quantity 

(pi + ••• + p„) — 

is maximal for the equilateral triangle. When n = 1, every maximizer is 
equilateral. 

Analogous results hold among parallelograms and ellipses, with squares 
and disks being the maximizers, respectively. 

Schrodinger eigenvalues. Consider the Schrodinger eigenvalue problem 

-hAu + Wu = £u 

in the plane, with Planck constant h > and real-valued potential W G 
L'^^iM?) that tends to +oo as |x| —7- oo. The spectrum is discrete [49 [ 
Theorem XIII. 67], and the eigenvalues £j are characterized in the usual way 
by the Rayleigh quotient 

R[u] = ^^ — - — ' ue H\R^)nL\W). 

Jir2 dx 

Here L^(VF) denotes the weighted space with measure \W\ dx. 

Once more we show that a rotationally symmetric situation maximizes 
the sum of eigenvalues. 

Theorem 3.5. IfW has rotational symmetry of order greater than or equal 
to 3, then 

{£l-\ ^ ^n)\2h/l\T--^Wjjs,WoT~-^ < {£l-\ ^'^n)L,W 

for each n > 1 and each invertible linear transformation T of the plane. 

When n = 1, equality holds if and only if T is a scalar multiple of an 
orthogonal matrix. 

The subscript "/i, W" on the right side of the inequality specifies the 
potential to be used and the value of the Planck constant, and similarly for 
the subscript on the left side of the inequality. 

This Schodinger result formally implies the Dirichlet result in Theorem l3.1l 
by taking W = on D and W = +oo off D, and choosing h = ||r^^||^g/2. 
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More general quadrilaterals. Among quadrilaterals we have so far han- 
dled only the special class of parallelograms. Now we show how to handle a 
larger class of quadrilaterals having two "halves" of equal area. 

To construct such domains, first write the upper and lower halfplanes as 

= {{xi,X2) : X2 > 0}, Rl = {(xi,X2) : X2 < 0}. 

Choose two linear transformations T+ and T_ that agree on the xi-axis, 
with T± mapping onto itself. Then the map 



r(x) 



T+x ifx£Rl, 



defines a piecewise linear homeomorphism of the plane mapping the upper 
and lower halfplanes to themselves. Assume also detT+ = detT_, so that 
T distorts areas by the same factor in the upper and lower halfplanes. 

We will not need explicit formulas for the linear transformations and 
T_, but for the sake of concreteness we present them anyway: 



a c± 
b 



where a / 0, 6 > and c± G M. 

Let D be the square with vertices at (±1,0), (0, ±1). Our goal is to show 
that this square maximizes eigenvalue sums among quadrilaterals of the form 
E = T{D) . These quadrilaterals have two vertices on the xi-axis and have 
upper and lower halves of equal area. 

Write 

Io{E)= [ \x\'^dx 
Je 

for the moment of inertia about the origin, for a domain E. 



Theorem 3.6 (Quadrilaterals with equal-area halves). Let D be the square 
with vertices at (±1, 0), (0, ±1) . Then for every map T constructed as above, 



(Ai + --- + A„) — 

-'0 



^3 

< (Ai + ••• + A„) — 

T{D) 



D 



for each n > 1. 

The inequality holds also for Neumann eigenvalues. 

The moment /q is generally greater than the moment of inertia I for the 
domain T{D), because the centroid of T{D) need not be at the origin. For 
example, if both r+ and r_ are shear transformations towards the right, 
then the centroid of T{D) will lie on the positive xi-axis, to the right of 
the origin. Meanwhile, the centroid of the rotationally symmetric domain 
D will always lie at the origin, so that I{D) = Iq{D). 

We conjecture that Theorem 13.61 can be strengthened to use I instead of 
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For the first eigenvalue (n = 1), Freitas and Siudeja [l5] showed recently 
with a computer-assisted proof that 

^1 ^2 ^~ ^3 ^~ ^4 

is maximal for rectangles among all quadrilaterals, not just among quadri- 
laterals with halves of equal area. For parallelograms, this result and The- 
orem [321 give the same information (see formula (j3.6p ). For general quadri- 
laterals we cannot easily compare the two results, because the relationship 
between the sum of squares of side lengths and the moment of inertia is 
unclear. 

4. Open problems for general convex domains 

For the Dirichlet fundamental tone we raise: 
Conjecture 4.1. Suppose $7 is a bounded convex plane domain. Then 

< 127r^ 



.2 - A' 



with equality on the right for equilateral triangles and all rectangles, and 
asymptotic equality on the left for degenerate acute isosceles triangles and 
sectors. 

The convexity assumption is necessary on the right side of the conjecture 
because otherwise one could drive the eigenvalue to infinity without affecting 
the area or moment of inertia, by removing sets of measure zero (such as 
curves) from the domain. 

The maximizer cannot be the disk in the last conjecture because triangles 
and rectangles yield a larger value, as we observed already in (|3.5p . 

As evidence for the conjecture, we note that AiA^// is bounded above and 
below on convex domains by an Inclusion Lemma, as was shown by Polya 
and Szego \A8\ §1.19, 5. lib]. They further evaluated XiA"^/! for a variety 
of triangles, sectors, degenerate ellipses and degenerate sectors [ISl p. 267]. 
Asymptotic expansions can be obtained also for degenerate triangles |14j. 
We examine the family of isosceles triangles in Figure [2j All this evidence 
is consistent with Conjecture 14.11 

For the first nonzero Neumann eigenvalue, we know ^12-^ /I is definitely 
maximal for the disk among all bounded domains, by an inequality of Szego 
and Weinberger (see Section [8]) . This quantity has no minimizer because it 
approaches zero for a degenerate rectangle. 

Now consider sums of eigenvalues. 

Conjecture 4.2. Suppose Q is a bounded convex plane domain. Then for 
the Neumann eigenvalues, 

(^2 + ---+/in) — 
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127r2 




67r2 



97rV2 







7r/3 



7r/2 



Figure 2. Numerical plot of the normalized Dirichlet fun- 
damental tone \iA^/I for isosceles triangles of aperture 
a G (0,7r). The maximizer is equilateral (a = vr/3), and 
the minimizer is degenerate acute (a — )• 0). 

is maximal when Q is a disk, for each n >2. 

The conjecture is true for the special case of ellipses by Corollary 13.21 
For Dirichlet eigenvalues, the conjecture fails because the square gives a 

larger value than the disk for (Ai+- • •+A„)^'^// when n = 1, 2, 3, 5, 6, 9, 10, 12; 

the disk does give a larger value for all other n < 50, and we suspect for 

n > 50 as well. 

5. Consequences of symmetry: tight frames, and moment 



In this section we recall the tight frame property of rotationally symmetric 
systems of vectors, and develop a moment of inertia formula for the linear 
image of a rotationally symmetric domain. These elementary consequences 
of symmetry will be used in proving Theorem 13. 1[ 

Tight frames. Let > 3 and write Um for the matrix representing ro- 
tation by angle Inm/N, for m = 1,. . . ,N. For each nonzero y £ M?, the 
rotations generate a rotationally symmetric system {Uiy, . . . ,UNy} in the 
plane. For example, the system consists of the A^th roots of unity when 



We start with a well known Plancherel-type identity for such systems. 
Lemma 5.1. Let N >3. For all column vectors x,y £M? one has 



matrices 




1 



N 




N 



Y,\x-{Uray)\^ 



m=l 
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Proof of Lemma \5.1[ We may suppose x and y have length 1 and He at 
angles 9 and (p to the positive horizontal axis, respectively. Then 

N N 
- ^ |x • = ^ 5] cos2(0 - - 2vrm/iV) 

m=l r?i=l 

1 1 ^ 

= 2 + 2iV ^ """'^^^ " " 2WiV) 

m=l 

ni=l 

_ 1 

~ 2 

as desired. The assumption > 3 ensures that e"^^'^/^ ^ 1 when summing 
the geometric series, in the last step. □ 

Figure [3] illustrates the lemma for = 3: it shows the projection formula 
Ylm=i(^ ' Umy)Umy = §2^ for a typical x G R^, where y = ( ^ ) is the vertical 
unit vector and Um denotes rotation by 27rm/3. Dotting the projection 
formula with x yields the Parseval-type identity in Lemma l5.1i 




Figure 3. The "Mercedes-Benz" tight frame (A^ = 3) in the plane. 

The lemma says that the rotationally symmetric system {C/iy, . . . , Unv} 
forms a tight frame. Readers who want to learn about frames and their 
applications in Hilbert spaces may consult the monograph by Christensen 
[To) or the text by Han et al. [19]. 

We next deduce a tight frame identity in which the vector y is replaced 
by a matrix. 
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Lemma 5.2. Let N > 'i,K > \. For all row vectors x and all 2 x K 

real matrices Y one has 

N 



m=l 

Proof of Lemma \5. 21 Write yi, . . . , yx for the column vectors of Y, so that 
\xUmY\'^ = Ylk=i l^^mUkl'^- Now apply Lemma EU □ 

Hilbert Schmidt norms and moment of inertia. When proving Corol- 
lary [321 we will need to evaluate the Hilbert-Schmidt norm of T^^ in terms 
of moment of inertia and area. 

Lemma 5.3. // the bounded plane domain D has rotational symmetry of 
order N > 3, and T is an invertible 2x2 matrix, then 

Proof of Lemma 15. The centroid of D lies at the origin, in view of the 
rotational symmetry of D. Thus the centroid of TD also lies at the origin. 

The moment matrix of D is defined to be M{D) = [Jj^ XjXk dx]j^k- We 
show it equals a scalar multiple of the identity, as follows. Let U denote 
the matrix for rotation by 27r/N. The rotational invariance of D under U 
implies that M{D) = UM{D)U\ so that if X is an eigenvector of M{D) 
then so is Ux, with the same eigenvalue. Since x and Ux span (using 
here that > 3), we conclude every vector in is an eigenvector with 
that same eigenvalue. Thus M{D) is a multiple of the identity. 

In particular, the diagonal entries in M{D) are equal. Since they sum to 
the moment of inertia L{D), we have 

M(D) = i/(D)(j J). (5.1) 

The moment of inertia of T{D) can now be computed as 
L{TD) = trM{TD) 

= tTTM{D)T^detT\ 

= h{D) ( tr TT^) I det T\ by dO]) 

= lL{D)\\T\\ls\detT\. (5.2) 

This formula gives us the Hilbert-Schmidt norm of T, whereas we want the 
Hilbert-Schmidt norm of T~^. Fortunately, the two are related, with 

r"'llH5 = rilW|detr|2 (5.3) 

by the explicit formula for T^^ in terms of the matrix entries, in two dimen- 
sions. Hence 

I{TD) = ^I{D)\\T-YHs\detTf, (5.4) 
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from which the lemma follows easily. □ 

An interesting consequence of the last lemma is that the moment of inertia 
of a linear image of a rotationally symmetric domain equals the moment of 
inertia of its inverse image, after normalizing by the area. 

Lemma 5.4. // the bounded plane domain D has rotational symmetry of 
order N > 3, and T is an invertible 2x2 matrix, then 

Proof of Lemma [5^71 By ()5.2p . and then using (j5.4p with T replaced by T~^, 
we find 

I{TD) _ 1 I{D) \\T\\%s _ I{T-^D) 
A{TDY ~ 2 A(L>)2 |detr| ~ Kf^Wf' 

□ 

The lemma holds also with T^'^ D instead of T^^D, since T"^ and T have 
the same Hilbert-Schmidt norm and determinant. 

6. Proof of Theorem 13.11 

We prove the Dirichlet case of the theorem. The idea is to construct trial 
functions on the domain T{D) by linearly transplanting eigenfunctions of 
D, and then to average with respect to the rotations of D. The Neumann 
proof is identical, except using Neumann eigenfunctions. 

Let ui,U2,U3, . . . be orthonormal eigenfunctions on D corresponding to 
the Dirichlet eigenvalues Ai, A2, A3, . . .. Consider an orthogonal matrix U G 
0(2) that fixes D, so that U{D) = D. Define trial functions 

Vj = Uj o U o 

on the domain E = T{D), noting vj £ Hq{E) because Uj G Hq{D). 
The functions vj are pairwise orthogonal, since 

/ VjVkdx= / UjUk dx ■ \ detTU^^l = 
J E Jd 

when j 7^ k. Thus by the Rayleigh-Poincare principle, we have 

j=l j=l JE «^ 

For each function v = vj we evaluate the Rayleigh quotient as 
fj^\Vv\'^dx _ /p|(Vn)(x)C/r-i|2(ix- |detr;7-i| 
f^v'^dx fjju'^dx-\detTU~^\ 

\{Vu)UT~^\'^ dx, (6.2) 

D 

where the gradient Vu is regarded as a row vector. In the last line we used 
that u = Uj is normalized in L'^{D). 
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Since D has A^-fold rotational symmetry for some > 3, we may choose 
U to be the matrix Um representing rotation by angle 27rm/A^, for m = 
1, . . . , A^. By averaging (j6.ip and (j6.2p over these rotations we find 

j^HE) <Y,j {if; |(Vn,)f/„r-ip}dx 

j=l j=l ''^ m=l 

" /" 1 

= {-\Vuj\'^\\T-'^\\'jjs} dx byLemmaEa 

n 

= l^\\T''\?HsY.^3iD), 

i=i 

which proves the inequality in Theorem 13.11 

Equality statement for Dirichlet fundamental tone, n = 1. Suppose equality 
holds in the theorem for the first Dirichlet eigenvalue. That is, suppose 

X,{T{D))=^\\T-YHS^iiD)- (6.3) 

We reduce to T being diagonal, as follows. The singular value decompo- 
sition of T can be written T = QRS where Q and S are orthogonal matrices 
with det S = 1 (so that 5" is a rotation matrix) and R= ( ) is diagonal 
with ri,r2 > 0. If ri = r2 then T is a scalar multiple of an orthogonal 
matrix. So suppose from now on that ri 7^ r2. 

Write D = S{D), so that D has rotational symmetry of order N. Note 
\i{D) = Ai(5), and that T{D) = QR{D) so that 

Xi{T{D))=X,{R{D)). 

Also 

11^ """IIhS = 11-^ """IllfS = ■'^1 ^ + ^2 ^• 

Hence equality in (j6.3p implies 

Ai(i?(5)) = i(rr2 + r2-2)Ai(5), 

which means that equality holds in (j3.ip for the domain D under the diagonal 
linear transformation R. 

Write u = ui for a first Dirichlet eigenfunction on D, so that 

Mxixi + Wx2X2 = -Ai(-D)it. (6.4) 

Inspecting the proof of the theorem, above, we see that one of the trial 
functions on R{D) is v = uoR^^, in other words f (xi, X2) = ^(a^i/ri, ^2/^2). 
Since equality holds in the Rayleigh principle (j6.ip with n = 1, we deduce 
that this trial function must actually be a first eigenfunction on R(D). That 
is, 

Av = -\i{R{D))v, 



16 

which means 

r-f^n^ixi + r2^u-,2X2 = -^(^r^ + r2^)M{D)u. (6.5) 

By solving the simultaneous linear equations (j6.4p and ()6.5p (which is pos- 
sible since ri 7^ we find that 

Uxixi = Ux2X2 = (6.6) 

This last formula must apply also if we rotate u through angle 2tt/N, because 
that rotate of u was used in one of the trial functions in the proof of the 
theorem above. Hence the second directional derivative of u in direction 
9 = 2tt/N must equal —^Xi{D)u. That second derivative is 

(cos^ 6')iixixi + 2(cos sin 6)u.c^x2 + (sin^ (^)ux2X2 , 

which equals —^Xi{D)u+sm{29)uxiX2 by (|6.6p . We conclude sin(2^)u^^2,'2 = 
0. 

Suppose N ^ A. Then sin(26') = sin(47r/A^) 7^ 0, and so Uxix2 = in D. 
Then u = Fi{xi)+F2{x2) for some functions Fi and F2, and substituting this 
formula into ([ES]) gives that F{'{xi) = -iAi(5)[Fi(xi) + ^2(^2)]. Taking 
the X2 derivative shows that F2 is constant. Similarly Fi is constant, and 
so u is constant, an impossibility. 

Therefore = 4. Write uj = \J Xi{D)/2. The equations (j6.6p say Ux^^xi = 
Ux2X2 = — w^u, and so 

u(xi, X2) = A cos(a;xi) cos(a;x2) + B sin(a;xi) sin(a;x2) 

+ C cos(a;xi) sin(a;x2) + -D* sin(a;xi) cos{ujX2) 

on D, for some constants A,B,C,D^. The 4-fold rotational symmetry of 
the domain further implies that each of the four terms 

j4 cos(a;xi) cos(a;x2) (6.7) 

i? sin(a;xi) sin(a;x2) (6.8) 

(C^ + Dl) cos(a;xi) s].n{ujX2) (6.9) 

(C^ + Dl) sin(a;xi) cos{ujx2) (6.10) 

is by itself a Dirichlet mode for D with eigenvalue 2a;^ = Xi{D), or else 
is identically zero, as we will now show. First, by adding and subtracting 
ti(xi,X2) and n(— xi, — X2) (its rotation by vr) we find that the functions 

/(xi, X2) = Acos(a;xi) cos(a;x2) + B sin(a;xi) sin(a;x2) 

and 

g{xi,X2) = C cos(a;xi) sin(ci;x2) + sin(a;xi) cos(a;x2) 

are each eigenfunctions on D (or else are identically zero). By adding and 
subtracting /(xi, X2) and /(— X2, xi) (rotation by 7r/2) we find that (|6.7p and 
(|6.8p are each eigenfunctions (or else are identically zero). By considering 
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Cg{xi,X2) — D^:g{—X2,xi) and D^:g{xi, X2) + Cg{—X2, xi) we learn that ()6.9p 
and ()6.10p are each eigenfunctions (or else are identically zero). 

The fundamental Dirichlet mode does not change sign. The nodal do- 
mains for each of the functions (|6.7p - (|6.10p are squares, and so D must lie 
within one of those squares. For (j6.8p . ()6.9p and ()6.10p . rotation by angle 
vr maps each nodal square to a completely disjoint square, which means 
that D cannot have 2-fold rotational symmetry, let alone 4-fold symmetry. 
Hence (j6.8p . (|6.9p and (|6.10p must not be eigenfunctions, and so necessarily 
B = C = D^: = 0. Thus the eigenfunction is (|6.7p . Taking ^4 = 1, we have 

u = cos{ujxi) cos(a;x2). 

Rotation by vr rules out every nodal square except the one centered at the 
origin, which is (— 7r/2w, 7r/2a;)^. Hence D is contained in this square. 

The square has first Dirichlet eigenvalue 2lo'^, which equals Ai(L'). Thus 
D must fill the whole square (except perhaps omitting a set of capacity 
zero, which does not affect the fundamental tone [E]). Then R{D) is a 
rectangle, and D = S~^{D) is a square and T(D) = QR{D) is a rectangle. 
This completes the proof of the "only if" part of the proof of the equality 
statement. 

For the "if" part of the equality statement, suppose D is a square and 
T{D) is a rectangle (possibly with sets of capacity zero removed). By ro- 
tating and reflecting D and T[D) suitably, we can suppose they have sides 
parallel to the coordinate axes and that T = ( ) for some ri,r2 > 0. 
Writing L for the side length of the square, we have 

\i{D) = 2(vr/L)2, \i{T{D)) = {tt /nLf+{7: /r2Lf = ^||r-i||2,^Ai(I)), 
so that equality holds in (j3.ip with n = 1. 

Equality statement for second Neumann eigenvalue, n = 2. Suppose equal- 
ity holds in the theorem for the second Neumann eigenvalue. Most of the 
preceding argument in the Dirichlet equality case applies without change, 
simply replacing Ai with fj,2 and the Dirichlet eigenfunction ui with the 
Neumann eigenfunction U2, and replacing the word "Dirichlet" with "Neu- 
mann". The argument works because the first Neumann eigenvalue of D 
is zero, with constant eigenfunction ui = const., and so the trial function 
vi = ui oUoR^^ is also constant and hence is a first eigenfunction on R{D). 
Thus if equality holds in the Rayleigh principle (16. ip for n = 2 then the trial 
function V2 = U2 o U o R~^ is a second eigenfunction on R{D). 

The significant difference from the Dirichlet proof begins at the sentence 
"The fundamental Dirichlet mode does not change sign." The second Neu- 
mann eigenfunction u = U2 does change sign on D: it has exactly two nodal 
domains {u > 0} and {u < 0}, each of which is connected. (We know the 
eigenfunction has at least two nodal domains because u is orthogonal to 
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the constant eigenfunction; it has at most two by Courant's nodal domain 
theorem (Tj p. 112].) 

Consider each of the four possible forms of u in turn, namely (j6.7p - (j6.10p . 
Each one has square nodal domains, and the two nodal domains of u in 
D must be subsets of such squares. Hence D intersects exactly two of 
the squares. At the same time, D has 4-fold rotational symmetry. These 
requirements prevent (16. 7p from being an eigenfunction for D, because if D 
intersected two of the nodal squares, then it would have to intersect at least 
five of them. Hence ^ = 0. Similarly ()6.8p cannot be an eigenfunction, and 
so B = 0. 

Next we deal with (j6.9p . (The argument is similar for (j6.10p .) Suppose 
+ > in (|6.9p . so that we may take 

M = cos(a;xi) sin(a;x2). (6-11) 

Then in order for D to intersect exactly two of the nodal squares, they must 
be the squares adjacent to the origin, so that 

D C (-7r/2tJ,7r/2w) x (-tt/cj, tt/cj). (6.12) 

We will deduce a contradiction below, so that necessarily + = 0. 
Hence none of the functions (j6.7p ~ (j6.10p is an eigenfunction, and so the case 
A'^ = 4 cannot occur. Therefore the only way for equality to hold is to have 
fi = f2, so that T is a scalar multiple of an orthogonal matrix. 

To obtain the desired contradiction, we will examine how the Neumann 
boundary condition is affected by the linear transformation. Since the do- 
main D has Lipschitz boundary, there exists an outward normal vector 
(ni,n2) at almost every boundary point (with respect to arclength mea- 
sure). At each such point {xi,X2) G dD, we know u satisfies the Neumann 
(or natural) boundary condition 

= Vu • (ni,n2) = Uxj^ni + u^^m] 
here we used that u as defined by (16. lip is globally smooth. Further, the 
point (riXi,r2X2) G dR{D) has an outward normal vector (ni/ri, n2/r2). 
Since v{xi^X2) = u{xi/ri^X2/r2) is a smooth eigenfunction on the closure 
of R{D), it satisfies the Neumann boundary condition: 

= Vv- (ni/ri,n2/r2) = Ux^ni/rl + Ux2n2/rj. 

Recalling that ri / r2, these simultaneous equations imply 

UxiiT-i = and Ux2'n2 = 0. 

Hence either Ux^ = or = 0, a.e. on dD. Recalling the formula (j6.1ip 
for u and the constraint ()6.12p on Z), we deduce that almost every boundary 
point is contained in the lines {xi = 0, ib7r/2a;}, {3:2 = 0, ib7r/2a;, itvr/a;}. 
Furthermore, we can rule out the vertical lines {xi = ib7r/2a;} because on 
those lines Ux^ 7^ and so ni = 0, which means the normal would be vertical 
and the tangent horizontal, so that the boundary would depart the given 
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vertical lines. Similarly we rule out the horizontal lines {x2 = Ojin/uj}. 
Hence the boundary of D must lie in the union of the lines {xi = 0}, {x2 = 
±7r/2a;}. Since these lines fail to bound a domain, we have arrived at a 
contradiction, as desired. 

Why did Polya not prove our theorem? Polya proved Theorem 13. II for 
the first Dirichlet eigenvalue Ai, except that he proved no equality statement. 
His result appeared in [44| and its proof in [47', Chapter IV]. 

Why did he not prove the theorem for sums of eigenvalues, or for Neumann 
eigenvalues, as we do in this paper? Or for higher dimensions as we do in a 
forthcoming paper [35]? 

A possible reason is that our method is subtly different from Polya's. 
We use rotational symmetry at a later stage in the argument. This delay 
permits us to handle more than just the first eigenvalue, and to handle 
Neumann eigenvalues too. Let us explain in more detail. Polya began by 
using rotational symmetry of the domain to obtain rotational symmetry 
of the fundamental Dirichlet eigenfunction ui: he observed that the rotate 
of ui is itself a positive eigenfunction and so must equal ni. Then Polya 
deduced that 

JoydxiJ JD\dx2j 2Jf) JDdxidx2 

Hence his linearly transplanted trial function ui oT~^ has Rayleigh quotient 
/^|VKoT-^)pdx _ /^|(Vm)(x)r-ipdx _ l,,y-i,,2 /plV^iP^a; 

as desired. 

The difficulty when trying to extend Polya's approach to sums of eigen- 
values is that the higher eigenfunctions are usually not symmetric under 
rotations, because of sign changes. The insight that permits us to prove 
Theorem 13. II is that while the rotate of a higher eigenfunction need not equal 
itself, it must still be an eigenfunction with the same eigenvalue, and thus 
can still be used to generate trial functions by linear transplantation. Our 
proof uses the whole family of rotations to generate many trial functions, 
and then averages over the resulting family of inequalities. This approach 
applies (without change!) to the Neumann eigenvalues too. 

7. Proofs of other results 

Proof of Corollary 13. 2L Every triangle can be written (after translation) 
as the image under a linear transformation T of an equilateral triangle cen- 
tered at the origin . Hence the inequality for triangles in Corollary 13. 2l follows 
from Theorem 13.11 and Lemma 15.31 The statements about parallelograms 
and ellipses are proved similarly. 
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Remark on Dirichlet maximizers when n > 2. It is not clear how to de- 
termine all maximizing domains for sums of eigenvalues beyond the first. 
For example, some (but not all) non-square rectangles can maximize (Ai + 
■ ■ ■ + ^n)^^/I when n > 2, as we now show. Consider a rectangle with 
side lengths li,l2, so that the area is A = I1I2, the moment of inertia is 
I = hl2{ll + ll) / 12 and 

_ 12 

T~ i-^ + i-^- 

The fundamental tone is 

Ai = 7r2(Z-2 + /-2). 

Notice Aivl^// = 127r2 for every rectangle (not just for the square), as we 
have observed before. Thus every rectangle is a maximizer when n = 1. 

Now fix n > 2, and fix the side length I2. For li sufficiently large, we have 
eigenvalues Xj = 7T'^{j'^li'^ + Z^^) for j = 1, . . . , n, and so 

lim (Ai + • • • + A„)— = 127r2n. 
Meanwhile, the square satisfies 

A^ A^ 
(Ai H h A„)— > nAi— = 127r^n. 

Hence for sufficiently large li, the rectangle with side lengths li and I2 is not 
a maximizer. 

Nonetheless, the rectangle can be a maximizer for some values of h and 
n. For example, let n = 3 and suppose the side lengths of the rectangle 
satisfy I2 < h < a/8/3/2- Then by simple comparisons we find 

Ai = 7r2(/-2 + Z2-2), 

A2 = 7r2(2V + «2"')> 
A3 = 7r2(Z-2 + 2%-2), 

and so 

(Ai + A2 + A3)^ = 727r2. 

This value is the same as achieved for the square {h = I2), and so there are 
many non-square maximizers when n = 3. 

The idea behind this construction is to identify a range of (/i,/2) values 
for which the eigenvalues Ai, A2, A3 have values Tr'^{j'^l^'^ + k'^l2'^) for (j, k) = 
(1, 1), (2, 1), (1, 2). This set of index pairs in 7? is invariant with respect 
to interchanging j and k. Hence Ai -|- A2 + A3 is proportional to +I2 , 
which allows us to cancel the denominator in A? /I and obtain an expression 
independent of the side lengths. This construction can of course be extended 
to arbitrarily large values of n, if desired. 
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Proof of Theorem 13. 3L The proof goes exactly as for the Dirichlet and 
Neumann cases in the proof of Theorem 13.11 except that for the Robin 
eigenvalues we must take account also of a boundary integral in the Rayleigh 
quotient: 

faEV^dsix) ^ jQ^u{UT-'xfds{x) 
J^v'^dx J^u{UT~^xy dx 

_ Jg^u{Ux)''\TT{x)\ds{x) 
u{U x)2 dx ■ I det T| 

by X I— )• Tx, where t(x) denotes the unit tangent vector to dD at x. Geo- 
metrically, |rr(x)| is the factor by which T stretches the tangent direction 
to dD at X. 

The symmetry of D implies that the tangent vectors rotate according to 
t(U~^x) = U^^t(x), and so replacing x with U'^x in the last integral gives 

IaKf-^ = \detTr' [ uixf\TU'\ix)\ds{x). 

Choose U to be the matrix Um representing rotation by angle 27rm/N, 
for m = 1,...,N. Averaging the preceding quantity over m and apply- 
ing Cauchy-Schwarz gives the upper estimate 

\detT\-' l^^u{x)'[^ ^ \TU~Mx)\'y' ds{x) (7.1) 



m=l 



= |detr| "^^p||r||//5 / u{x)'^ ds{x) 
V2 JdD 

by Lemma [5T2| since |t(x)| = 1. Multiplying by a\\T'^\\Hs/V^ gives 

hlT'^Wns'^ I u{xfds{x) 

^ JdD 

by dO]). 

With the aid of this last estimate we can straightforwardly adapt the 
proof of Theorem 13. II to the Robin situation, and then call on Lemma [5.3l to 
interpret the Hilbert-Schmidt norm of in terms of moment of inertia. 

Equality statement for Robin fundamental tone, n = 1. The proof of the 
equality statement follows the Dirichlet case in Theorem 13. II up to the point 
where = 4 and 

u = cos(cjxi) cos(a;x2), 

and D contained in the open square S = (— 7r/2a;, 7r/2w)^. We want to 
deduce a contradiction, so that the only way for equality to hold when 
n = 1 is for T to be a scalar multiple of an orthogonal matrix. 

Equality must hold in the application of Cauchy-Schwarz at (j7.ip . except 
using R instead of T and using D instead of D. Hence 

|i?C/f V(x)| = \RU^^t{x)\ = |i2C/3"V(x)| = \RU^'^t{x)\ (7.2) 
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for almost every (with respect to arclength measure) x £ S dD; here we 
use that u{x)'^ > on 5. 

Consider such an x-value and write ri and T2 for the components of the 
tangent vector t{x). Then \R{ll)\ = \R {'^l^) \ hy ([72]), or 

{nnf + (rara)^ = (-rirs)^ + (rari)^. 

Since rf 7^ r^, we can simplify to = tI • Thus the tangent line at x has 
slope ±1, and hence so does the normal vector. 

The four possible normal vectors are n{x) = {ei, 82)/ where £1,62 £ 
{ — 1, 1}. Thus the Robin boundary condition ^ + au = says 

eiiia:i + e2Ux2 + V2au = 0. 
Substituting u = cos(a;xi) cos(cjX2) yields that 

£1 tan{ujxi) + £2 tan(a;x2) = V^a/u. 

We conclude that every point x €z S D dD lies on one of these four curves. 

These curves have slope ±1 at only finitely many points in the square S, 
and so we conclude that no points of dD lie in that square. Hence dD lies 
entirely in the boundary of the square S. The Robin condition fails on dS, 
though, because u = there while |^ / 0. 

This contradiction completes the proof. 

Proof of Corollary 13.41 First we show that among all triangles G, the 
quantity 



(pi H \- Pn 



(7.3) 

G 



is maximal for the equilateral triangle, for each n > 1. For this, let D be an 
equilateral triangle centered at the origin, and let G = T[D) where T is an 
invertible linear transformation. Note that by Lemma l5.3| 



^ / I{D) \\T-^\\hs ^ ^ I I{TD) _ ^ / I{G) 



A{Df V2 ^ A{TDf \ A{Gf 

Hence replacing a in Theorem O with a.JI{D)/A{Df proves that expres- 
sion (17. Sp is maximal when T is the identity, that is, when G is equilateral. 

Next, assume the triangle G has the same area as the equilateral D. Then 
the equilateral has smaller moment of inertia, I{D) < I{G), as can be proved 
from formula (jS.Sp for the moment of inertia along with Cauchy-Schwarz 
and the isoperimetric theorem for triangles. Hence by (j7.3p the expression 

A^ 

is maximal when G equals the equilateral triangle D. (Here we used mono- 
tonicity of the Robin eigenvalues with respect to the Robin parameter.) 
Finally, replacing a with a\J A(D)^/I{D) proves the corollary, for triangles. 



(pi H \- Pn 
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Argue similarly for parallelograms and ellipses. 

Proof of Theorem 13.51 The proof proceeds as for the Dirichlet case in 
Theorem 13. 1|. except that we must consider also the potential term in the 
numerator of the Rayleigh quotient. The key estimate is that 

f^{WoT-^)v^dx _ fjj{W o U-^)u^ dx ■ \detTU-^\ 

dx u'^ dx • \ det TU~^\ 



Wu^ dx 



D 



by the rotational symmetry of W . The proof is now easily completed. 

Incidentally, the assumption in the theorem that the potential W should 
grow at infinity can be significantly weakened [36j . 



Equality statement for Schrddinger fundamental tone, n = 1. Just like in 
the Dirichlet case, the singular value decomposition allows us to reduce to 
T being diagonal. The analogues of equations (j6.4p and (j6.5p are that 



Hux^xi + UX.2X2) = {W - £i{W))u, 

2h 



-2 , -2 
"1 +^2 



[ri^u^^xi + r^^u^^x^) = {W- £i{W)) 



u. 



(These equations hold pointwise a.e. by elliptic regularity theory, since the 
potential is locally bounded [171 Theorem 8.8].) Solving these simultaneous 
equations, we deduce (since ri / r2) that 

1 



^XiX\ — '^X2X2 



{W-£i{W))u. (7.4) 



2h 



The potential W{x) is assumed to tend to 00 as \x\ — )• 00, and so W — £i > 
whenever \x\ is sufficiently large. Multiplying (j7.4p by u and integrating in 
the xi direction, we deduce that — J^Uxi dxi > when \x2\ is sufficiently 
large, so that u{xi,X2) = for almost every xi. Since (I7.4p says that u 
satisfies the one dimensional wave equation with X2 playing the role of time 
variable and xi playing the role of space variable, we conclude that u = 
a.e. in M^. This contradiction completes the proof. 

Proof of Theorem 13.61 We prove a generalization of Theorem 13. 11 namely 
that 

(Ai + • • • + Xr,){TD) < \{\\T-Yhs + WTiYhs) i^i + ■■■ + ^n){D) (7.5) 

for any bounded D having rotational symmetry of order > 4 with N even. 

The proof of Theorem 13.11 requires some modifications. First we show 
that pairwise orthogonality of the vj remains valid. Decomposing D and 
E = T{D) into their upper and lower halves D± = DnM^ and E± = Ef^R\, 
we compute 

/ VjVk dx = I UjUk dx ■ I det T±U~^\ . 
Je± Jd± 
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These upper and lower terms sum to zero because det T+ = det T_ and 
Jj-, UjUk dx = by assumption, when j ^ k. Thus vjVk dx = 0. 
Next we consider the Rayleigh quotient of v. We decompose it as 



\Vv\^dx 



Y,! \{Vu){x)UT^^\^ dx, (7.6) 

^ JU(D+) 



f^v^dx Ju(D±) 

where in this calculation we use once more that the determinants of Tj^ and 
T_ agree. 

Since is even, C/7V/2 represents rotation by vr, so that Um+Nj2[D±) = 
Um[D^) and ?7m+Af/2 = —Um- Hence when we average (|7.6p over the rota- 
tions U = Um we obtain 

N 

.±^?dx 



\{yu){x)UmT^ 
N/2 

= T.m12{ +/ )\iVu)ix)UmT^'\'dx 

N/2 

= -TtTI l(Vtx)(x)C/„T±i|2dx since C/„(Z?) = D 

N 

± •'^ rn=l 

by Lemma 15.21 Now complete the proof of (|7.5p by recalling u = uj and 
summing over j. 

Then the theorem follows from (j7.5p and the evaluation of the Hilbert- 
Schmidt norms in the next lemma. 

Lemma 7.1. Let T be the piecewise linear homeomorphism in Theorem VJ.bX 
If the hounded plane domain D has rotational symmetry of order N > A, 
with N even, then 

\mrHs+\\TirHs) = f,{TD)/^,{D). 

Recall Iq denotes the moment of inertia about the origin. 



Proof of Lemma 7.1. The moment integrals over the upper and lower halves 
of D agree, with 



/ XjXk dx = 
Jd+ Jl 



because Z)+ maps to -D_ under rotation by it (that is, x 1— t- —x). Here we 
use evenness of the order of rotation. 
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Hence the moment matrices satisfy M(Z)+) = M{D^) = M{D)/2. Since 
M{D) = \h{D) ( Q 5 ) , as shown in the proof of Lemma 15.31 we deduce 

MP±) = i/o(Z)) (J J) (7.7) 

Now the moment of inertia of TD about the origin can be computed as 
Io{TD) = trM{TD) 

= trM(r+L>+) +trM(r_L»_) 

= J^(trr±M(D±)rl) -IdetTil 

= llo{D) ^ ( tr T±tI) ■ I det r± | hy ^ 

± 

= ^Io{D) {\\T+fHs + WT-Wls) I detTil, 

where in the last step we used that det T+ = det T_ . 

The Hilbert-Schmidt norm of the inverse T^^ is related to the Hilbert- 
Schmidt norm of T± by (|5.3p , and so 

loiTD) = ^Io{D){\\T-YHS + \\TlYHs)\detT^f, 
from which the lemma follows. □ 

8. Literature on maximizing low eigenvalues under area, 
perimeter, in-radius or conformal mapping normalization 

This paper gives sharp upper bounds on the sum of the first n > 1 eigen- 
values, normalized by /I. To help put these results in context, we now de- 
scribe results and conjectures that apply to the low eigenvalues (n = 1, 2, 3). 

Dirichlet eigenvalues. The quantity XiA^/L"^ (where L is the perimeter) 
is maximal among triangles for the equilateral triangle, by work of Siudeja 
|50j . This result is stronger than Polya's upper bound ([33]) on Ai^V^; 
because AL"^/! = 36{h + I2 + hf/{ll + + if) by ([33]) and this ratio is 
maximal for the equilateral triangle (when = /2 = ^3)- 

Further, the normalized spectral gap (A2 — \i)A^ /L/^ is maximal among 
triangles for the equilateral, by more recent work of Siudeja [51| . and thus 
/L"^ is maximal for the equilateral also. Hence (A2 — XijJ^ /I and 
A2^'^// are maximal for the equilateral, which improves on Corollary 13.21 for 
n = 2. 

Among general convex domains, Ai^^/L^ is maximal for degenerate rect- 
angles by work of Polya |45] . That result differs from our Conj ecture 14 . 1 1 on 
Aij4^//, where the equilateral triangle should also be a maximizer. 

Turning now to the indradius R, it is easy to see for triangles (or any 
polygon with an inscribed circle) that A/L is proportional to R. Hence 
the preceding upper bounds for eigenvalues of triangles using can 
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be restated using a normalizing factor of inradius squared. In particular, 
XiR^ is maximal for the equilateral triangle. A more general result is due 
to Solynin [52J: among all A^-gons with an inscribed circle, XiR^ is maximal 
for the regular A^-gon. Of course, among general domains the maximizer of 
Ai-R^ is simply the disk, by domain monotonicity. 

For area normalization, Antunes and Freitas [2] Conj. 6.1] conjecture 
that the Faber-Krahn lower bound on XiA has a sharp upper analogue 
that includes an isoperimetric correction term: they conjecture that among 
simply connected plane domains, 



with equality for the disk and (in a limiting sense) for degenerate rectangles. 

Under a conformal mapping normalization, Polya and Schiffer proved 
lower bounds on sums of reciprocal eigenvalues 1/Ai + • • • + 1/ , with the 
disk being extremal [U]. Extensions to surfaces with bounded curvature 
were proved by Bandle [TJ p. 120], and to spectral zeta functions and doubly 
connected surfaces by Laugesen and Morpurgo [30] |3T] . 

Lastly, the scale invariant ratio A2/A1 is maximal for the equilateral tri- 
angle among acute triangles, by work of Siudeja [51]. The conjecture re- 
mains open for obtuse triangles. For general domains, this Payne-Polya- 
Weinberger functional is known to be maximal for the disk, by Ashbaugh 
and Benguria [5j. 

Neumann eigenvalues. Stronger inequalities are known than the one we 
found for ij,2A'^/I in (|3.7p (which is the case n = 2 of Corollarv 13. 2p . In 
fact, /i2^ is maximal for the equilateral triangle among triangles, and for 
the square among parallelograms, and for the disk among all bounded plane 
domains. The first of these stronger inequalities was proved recently by the 
authors [32^ Theorem 3.1]. The second, for parallelograms, is unpublished 
work of the authors. The third inequality is a result of Szego and Wein- 
berger [541 I55| . These inequalities for 122 A are stronger because A"^/! is 
maximal for the equilateral triangle among triangles, for the square among 
parallelograms, and for the disk among all domains. 

Our inequalities in Corollary 13.21 hold for all n > 2. In contrast, the 
stronger inequalities fail to extend to n = 3. The maximizing domains are 
instead somewhat elongated: the "arithmetic mean" (^2 + IJ'3)A seems to be 
maximal among isosceles triangles for an aperture slightly greater than vr/G 
(according to numerical work), rather than for the equilateral triangle with 
aperture 7r/3; and {^2 + ^3)^ seems to be maximal among parallelograms 
for the 2:1 rectangle rather than the square (see [6i §5] for comments on 
rectangles). The maximizer among convex domains is apparently not known. 
The only positive result is that the disk is maximal among 4-fold symmetric 
domains [6l §4]. Incidentally, it is open to maximize the geometric mean 
y/fj,2tJ'3A. The disk is conjectured to be extremal, by I. Polterovich. 
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Among convex plane domains, it is open to maximize H2L'^ ■ The disk is 
not the maximizer, because the equilateral triangle and the square have a 
larger value (in fact, the same value). The maximizer for /U2-C'^, where D is 
diameter, is known to be the degenerate obtuse isosceles triangle by work of 
Cheng P, Theorem 2.1], [33^ Proposition 3.6]. For the problems mentioned 
above, and for related conjectures on triangles, see [32l §IX]. 

Sums of reciprocal Neumann eigenvalues were minimized by Dittmar [TT] , 
under conformal mapping normalization. 

Lower bounds. Sharp lower bounds on Dirichlet eigenvalue sums for tri- 
angles are proved in a companion paper [34j, under diameter normalization. 
Lower bounds for the Neumann eigenvalue //2 are found in an earlier work 
|33j . References to other lower bounds can be found in those papers. 
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Appendix A. Eigenvalues of equilateral triangles, rectangles, disks 

The Dirichlet eigenfunctions of equilateral triangles were derived about 
150 years ago by Lame |281 pp. 131-135]. (See the treatment in the text 
of Mathews and Walker [39', pp. 237-239] or in the paper by Pinsky [43] . 
Note also the recent exposition by McCartin [40j.) Dirichlet eigenfunctions 
of rectangles and disks are well known too [7j. The eigenvalues are: 

{(167r^/9) [j^ + jii2 + iVl '■ Jiii2 ^ l} foi" an equilateral triangle of side 1, 

{tt^ [{ji/hf + {j2/hf] ■■ ji,i2 > 1} for a rectangle of side lengths li,l2, 

{jrn,p '■ m > 0,p > for the unit disk, 

where jm,p is the pth zero of the Bessel function . 
The Neumann eigenvalues are: 

{ (167r^/9) [j^ + Jii2 +^1] • Ji)i2 > O} for an equilateral triangle of side 1, 
{^^ [Ui/hf + {j2/hf] ■■ ji,i2 > 0} for a rectangle of side lengths /i, /2, 
{(jm,p)^ : m > 0,p > l} for the unit disk, 

where j'^ p is the pth zero of the Bessel derivative . See [3 HI] . 

The Robin eigenvalues of rectangles and disks can be found by separation 
of variables. The eigenvalues are known also for the equilateral triangle |42j . 
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